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Abstract 

We determine the size effect on the iattice thermai conductivity of nanoscale wire and multilayer 
structures formed in and by some typical semiconductor materials, using the Boltzmann transport 
equation and focusing on the Knudsen flow effect. For both types of nanostructured systems we 
find that the phonon transport is reduced significantly below the bulk value by boundary scattering 
off interface defects and/or interface modes. The Knudsen flow effects are important for almost all 
types of semiconductor nanostructures but we find them most pronounced in Si and SiC systems 
due to the very large phonon mean- free paths. We apply and test our wire thermal-transport 
results to recent measurements on Si nanowires. We further investigate and predict size effects in 
typical multilayered SiC nanostructures, for example, a doped-SiC/SiC/Si02 layered structure that 
could define the transport channel in a nanosize transistor. Here the phonon-interface scattering 
produces a heterostructure thermal conductivity smaller than what is predicted in a traditional 
heat-transport calculation, suggesting a breakdown of the traditional Fourier analysis even at room 
temperatures. Finally, we show that the effective thermal transport in a SiC/Si02 heterostructure 
is sensitive to the oxide depth and could thus be used as an in-situ probe of the SiC oxidation 
progress. 

PACS numbers: 44.05. +e, 44.10. +i, 66.70.+f 
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I. INTRODUCTION 



Thermal transport in nanostructured semiconductor systems is important for the contin- 
ued development of nanoscale electronics.—^ A number of recent experimental 3 - ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 and 
theoretical 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ' 17 !^ studies contributes to an 

emerging understanding of interesting science issues and potential technology problems. 
Opto- and power-electronics devices generate significant heat that must be dissipated to 
avoid degradation of the device performance,- 2 ^ and similar problems increasingly affect 
also more standard electronics applications. 

The thermal transport in semiconductor structures is dominated by phonons and mea- 
surements on bulk samples would indicate excellent thermal transport properties due to 
the material hardness and purity. However, typical device designs involve necessary nano- 
structuring, for example, the construction of material heterostructures that define quantum 
wells or traps a transport channel below an oxide and above a doped substrate. The associ- 
ated formation of interfaces is documented to effectively suppress the thermal transport both 

p ara V| c |3.4.5.6.10.11. 12.13.14.15.16.17.18.19.20.21.22. 23.24.25.26 an( j p cl -p cni jj cu l a r 9.27.28.29.30.31. 32.33.34.35 ^ Q 

the materials boundaries due to the phonon Knudsen flo w 10 ! 11 ' 15 ' 37 ' 38 and the phonon re- 
flectio n 6 ! 12 ! 28 ! 29 ! 34 effects, respectively. Both types of interface mechanisms lie outside the 
traditional Fourier analysis 3 ^ of heat conduction but produce a suppression of the effective 
thermal conductivity within a region extending to about one phonon mean-free path £ m f p 
from the interface. The interface impact of scattering can even cause crossovers in the 
nanostructure conductivity such that a material with high bulk thermal conductivity is not 
necessarily the best thermal conductor upon nanostructuring.— While thermal-conductivity 
suppressions can work to improve the efficiency of both thermoelectric and thermoionic 
cooling , 1 ! 40 ! 41 ! 42 ! 43 ! 44 ! 45 ! 46 room-temperature estimates,— £ m f p (T = 300K) > 100 nm, also tes- 
tifies to some adverse effects even for present devices. Worse, the quest to nanosize modern 
electronics can only exacerbate the thermal dissipation problems^ 6 - since the increase in 
packing density simultaneously will produce a further degradation in the effective thermal- 
transport properties. 

In this paper we illustrate the general nature of this Knudsen effect on the effec- 
tive lattice thermal conductivity of semiconductor devices by documenting the effects on 
two simple sample structures: thin semiconducting wir o 3 i 15 i 16 i 17 i 19 i 20 i 21 i 48 and multilayered 
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systems,— i 5 i 6 i 7 i 10 i n i 12 i 13 i 14 i 35 i 47 Fig. ^ Semiconductor wires and layers represent examples of 
semiconductor nanostructures that are interesting in themselves. They can also be viewed 
as typical building blocks in general semiconducting devices that typically have transport 
channels sandwiched between an oxide and substrate. 

Figure Q shows schematics of both types of nanostructures; the light regions in the 
schematics represent the transport active systems of high thermal conductivity while the 
dark regions illustrate the processed surrounding which may be oxides or heavily doped ma- 
terials. In both systems the oxide or the dopants will inhibit or limit the thermal transport 
in the processed regions and move the phonon transport to the transport active wire or 
layer. However, the interface scattering suppresses the thermal transport even in the (pure) 
transport channel where the scattering causes a Knudsen flow . 37 ' 38 

We obtain analytical results for the thermal-transport suppression given by special func- 
tions in both the wire and (multi)layer cases. The analytical results hold exactly when the 
phonon mean-free path £ m f p can be approximated as constant^ and simplify the thermal- 
transport calculations for wires and multilayers to simple frequency integrals in the more 
general case when £ m f p depends on mode Y and momentum 'g'. We thus extend previous 
investigations^ ll0lllll2ll3ll5ll8i:L9135 by efficient evaluations of the wire and layer thermal con- 
ductivities K wire , Ki ayer , and by (analytical) results for the layer Knudsen flow that arises 
when the nature of interface scattering differs at the two materials boundaries in a multi- 
layer structure. These results also permit an efficient evaluation of the multilayer thermal 
transport. We find that the Knudsen flow effect is important in nanosized wires and multi- 
layered structures formed from typical semiconductor materials due to the very large phonon 
mean-free paths. We test the accuracy of the approximations used to calculate the phonon 
Knudsen flow by comparing our results to experimental measurements for a nanosized Si 
wire.— We further investigate the phonon Knudsen flow in typical multilayer structures, in- 
cluding a simplified model of a nanosize transistor transport channel,^ 7 - and demonstrate a 
breakdown of the traditional Fourier analysis 3 ^ for some semiconductor nanostructures at 
and below room temperature. Finally we suggest an indirect but in-situ probe^ 7 - for the 
oxidation progress in SiC (and Si) materials. 

Figure El summarizes the predicted suppression in the lattice thermal conductivity of 
general semiconductor wire and layer nanostructures. We stress that there is a difference 
in thermal-conductivity range used to present the results for SiC, Si, and GaAs materials 
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but that the overall physical behavior and consequences for semiconductor nanostructures 
remain unchanged in all materials. The left panels show K wire (dashed curves) and «i a yer (solid 
curves) for varying thickness rf wirej i ayer for SiC, Si, and GaAs at T = 300 K. The dotted lines 
identify the experimentally observed value of the bulk-material thermal conductivities from 
which we also extract the listed values of the phonon mean-free paths £ m f P - The top right 
and middle right panels document the temperature variation in the Knudsen-flow thermal 
conductivity in SiC and Si layers (solid curves) and wires (dashed curves). Finally, the 
lower right panel shows the calculated temperature variation in the phonon mean-free path 
£mtp{T) for SiC, Si, and GaAs. The phonon Knudsen-flow effect on the thermal transport 
arises because this £ m f p (T) approaches the microscale at lower temperatures and may even at 
room temperature exceed the typical separation between boundaries in modern and future 
semiconductor devices. 

This paper is organized as follows. In section II we estimate the bulk-phonon mean- 
free path in SiC, Si, and GaAs, and we discuss the general Knudsen-flow behavior in low- 
dimensional systems. In section III we analyze the phonon-boundary scattering in wires and 
layered structures and we present complete special-function evaluation a 47l5 ° l51 for the lattice 
thermal conductivities K wire and /ti ay er- In section IV we discuss and compare the present 
approach with traditional thermal transport calculations of low- dimensional systems, and 
we suggest a possible in-situ probe of the SiC oxidation process. Finally, section V contains 
our conclusions. 



II. PHONON TRANSPORT THEORY 



We investigate the thermal transport properties in semiconductors where the dominant 
heat carriers are the phonons. For every phonon momentum q (and phonon mode i) we 
determine the distribution function N q ^ which relaxes to an equilibrium value, N®,^ = 
l/[exp(hujq /ksT) — 1], given exclusively by the temperature T and the phonon energy 
huif 1 of said mode. We use a relaxation time to describe the effective scattering in 
the bulk material. Further, we ignore the optical modes and assume for the longitudinal 
and transverse mode o 52l5S (i = l,t) a simple isotropic phonon dispersion given by q = 
uj(1 + aiU 3 )/ci and parameters ol{ fitted to the phonon energies at the zone boundaries 0j!*» 
Table I lists the sound velocities, Cj, and effective Debye temperatures, 0j, for SiC, Si, GaAs, 
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and Ge along with the fitted values for the dispersion parameter on. 

We mostly focus on the case with a constant, effective phonon mean-free path £ m f p (that 
only depends on temperature and materials) but we also provide results and analysis for the 
more general case of a (mode- and) momentum-dependent mean-free path £^ p (q) = v\ ,tq . 
We use our model of the phonon dynamics to estimate the group velocities Vq = q/(1 + 
4a ju; 3 ). In the following derivation and discussion we generally suppress the explicit labeling 
by mode index but all calculations have been performed by handling the longitudinal and 
transverse phonon modes separately. 

The distribution change SN q = N q — N2 induced by a small thermal gradient, VT, is cal- 
culated by solving the Boltzmann transport equation in the relaxation-time approximation,— 

(v, ■ V)N q = (1) 

to linear order in the thermal gradient VT. The q (mode) dependent value of the relaxation 
rate l/r g is effectively fitted from bulk-thermal transport measurements^*^ as discussed 
below. In turn, the calculated distribution change 5N q permits our evaluation of the thermal 
current 

j = £ / (0tK v ;^ (2) 

and hence of the thermal conductivity 

k = -i(VT)- 1 . (3) 

The interface scattering enters from the proper inclusion of the boundary conditions in 
Eq. (HJ). 



II A. Estimate of the phonon mean-free path in bulk 

To obtain an approximate evaluation of the phonon mean-free path we first consider 
the thermal transport in a bulk material where the solution of the linearized Boltzmann 
transport equation (1-BTE) simplifies to 

SNf^ = -rf (v« ■ VT) . (4) 

The lattice thermal conductivity, K Bu ik, is calculated by inserting this solution in Eqs. (J2J 
and 0. For a frequency and polarization-independent phonon mean-free path, Tq Vq = 

5 



£mip(T), the lattice thermal conductivity in a bulk material /tBuik oc £ m { p is simply pro- 
portional to this effective mean-free path . 11 ' 15 Comparison to experimental values for the 
temperature variation in the bulk thermal conductivities, Kg X ^ lk (T), therefore yields a simple 
phonon mean-free path estimat o 11 ' 15 

w^=^cn^» + ^r. (5) 

B \ t I / 

The estimate depends on the modal dispersion as expressed through the occupation- weighted 
integrals 

/.«. it 
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The lower right panel of figure 121 shows the calculated variation of the phonon mean- 
free path £ m{p as a function of the temperature for the set of semiconducting materials 



studied. For GaAs and Ge we use values of Kg X ^ lk taken from Ref. |52| while for SiC and 



Si we use values from Ref. The corresponding curve for Ge coincides almost exactly 
with that for GaAs because of the similarity in sound velocities and in measured values for 
the bulk thermal conductivities. We note that the phonon mean-free paths in the three 
materials are almost of micron-scale dimensions and that SiC possesses the largest value of 

f IfiSiC c«Si q /)Ge,GaAs\ 

Our approximative constant-£ m f p description provides simple and analytical evaluations 
of the Knudsen thermal-transport suppression in both the wire and multilayer cases. This 
description permits qualitative estimates of the general effect.^ The actual bulk-phonon 
mean-free path 4nf P (<?) = T q^ v (p wm > °f course, depend on the mode and on the phonon 
momentum q and/or the frequency uj q ^ and calls for a more refined analysis. We test the 
validity of the constant-£ m f p results by also considering the Knudsen flow arising with a 
frequency-dependent (mode-independent) effective decay rate^ 5 - r = A(T)/u 2 and by com- 
paring the results of both theory descriptions to recent experimental observations of the 
thermal transport in Si nanowires.— For the test model we determine the function A(T) 
from measurements 56 of the temperature variation of the bulk (Si) thermal conductivity 
using 

™ exD 67r 2 hf2g t (T) ^(T)\ _1 
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where the functions gi are defined as 




e t /r 



dx 



x A e x [l + ai{xk B T/h) 3 ] 2 
(e x -l) 2 [l + 4ai{xk B T/h) 3 Y 



(8) 



In this more refined model description we determine below the Knudsen suppression ana- 
lytically up to a single integration over the phonon frequency. 

II B. Phonon Knudsen flow 

We focus our study of the Knudsen flow effect on two types of nanosizes semiconductor 
structures that can be viewed as general building blocks in devices. The upper panel of 
Fig. 1 illustrates two possible realizations of wire systems^ l9ll5l4s which could be embedded 
in an oxide or could exists as a free-standing object. The phonons are assumed to undergo 
either specular or diffusive scattering with the probability p and 1 — p at the wire boundary. 
The lower panel of Fig. 1 shows a corresponding realization of a multilayer structure in 
which the transport- active channel (of high material purity) is located between two similar 
or different materials representing, for example, an oxide cab and a highly doped substrate. 
The surrounding materials (and/or vacuum) suppress the average transport not only because 
they have inferior conductivities themselves but also because their material boundaries cause 
a Knudsen-flow suppression even in the transport channel. Again, two types of scattering 
events are possible and phonons are assumed to scatter specularly from the upper (lower) 
interface with probability p+r_) and diffusively with probability 1 — respectively. The 

value of p + will generally differ from p_ as, for example, an oxide interface is rough (p + = 0) 
whereas the interface to an underlying (doped or undoped) substrate may be very smooth 



The phonon-boundary scattering off interface/surface defects becomes dominant when 



general situation we introduce q±_ (q\\) to denote the component of the phonon momentum 
q that is perpendicular to the material boundary Q (parallel to the thermal gradient, VT). 
Specular phonon scattering has only indirect effects on the overall transport because the 
phonons will have the same momentum component qu in the direction of the thermal gradient 
VT as before the scattering event. For such specular scattering we have 



(P 




the phonon mean-free path exceeds the dimensions of the material . 10 ! 11 ! 15 ' 38 To analyze the 



6N q (Q, q ± ) = 6N q (Q, -q ± ). 



(9) 
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In contrast, diffusive interface scattering causes a significant transport suppression as the 
phonon emerges with a random direction given by the local thermal distribution. The 
implicit thermalization ensures 

SN q (n, q ± >0) = 0, (10) 

and this boundary condition affects the distribution also away from the interface. 

To calculate the resulting impact on the thermal transport for our low-dimensional struc- 
tures we solve the 1-BTE 

fl 11 1 dN° 

V-V) + — \SN q = --— f (q v VT) (11) 

subject to the boundary conditions Eqs. Q and (fTU|) . We express the general solution 
relative to the corresponding bulk case 

5N q = 8Nf^[l-h q (r ± ,qJq,p)\ (12) 

using r_|_ to denote the coordinate direction away from the local surface. Furthermore, 
the form of the 1-BTE (II 1|) permits us to express the relative suppression in terms of an 
amplitude (fi q (q±/q,p) > and an exponential decay 

h q (r ± ,qjq,p) = <f )q (q ± /q,p)e~ r ^ v "^^ . (13) 

The amplitudes 4> q {q±/q,p) are fixed by analysis of the boundary conditions . 15 ! 37 ' 50 ' 51 

Figure El illustrates the general nature of the Knudsen flow in the phonon distribution, 
5N q /5Nl 3 ulk , as functions of the momentum angle qj_/q (top panel) and of the spatial coordi- 
nate x across the diameter (bottom panel). The Knudsen-flow distribution 5N q is specified 
by numerical integratio n 13 ' 15 ' 1 '? and/or analytical evaluation . 11 ' 47 ' 50 ' 51 The figure details the 
variation in 8N q / 5A^ ulk and compares results obtained assuming purely diffusive scattering 
(p = 0) for two wires with diameters d wire = £^j {p {q) and 2£^ (q). The calculated suppres- 
sion applies for all phonon momenta (and modes) for wires when, as we typically do, we 
approximate the phonon mean-free path as constant. More generally, the figure illustrates 
the universal effect in the Knudsen flow with a suppression of distribution changes extend- 
ing far away from the boundaries; similar momentum and spatial variation curves result 
both under more general assumptions for the mean-free path and in the case of multilayered 
semiconductor nanostructures. 
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The resulting effect on the thermal conductivity k can be summarized by averaging the 
relative suppression h q u\ over the cross section (perpendicular to VT). For every phonon 
q and i we define 7~L q i(q±/q,p) = J T dT h qi (r^,q^/q,p)/T and express the nanostructure 
conductivity change 

/j3 
— 4sKV - tN^H qi (qJq,p). (14) 
(2tt) q 

It is convenient to introduce an effective inverse Knudsen number r)i(q) = <i„i r e(iayer) / \ T l v q ) 
and describe the thermal-conductivity suppression^ 

6k . _ E» Jp' /T dxx 4 e x (e x - 1)~ 2 [l + ai(xk B T/h) 3 ] 2 F qt (p)\ vMx)) 



K Buik ^. J o Wl/J darxV (e x - l)" 2 [l + ai(xk B T/hf] 

as a weighted frequency integral over the effective suppression factor 

FMUi) = (y) J d\ (^j H gi (qjq,p). (16) 

We emphasize that significant simplification follows from an assumption of a mode and 
frequency independent £ m f p = constant. In that case we obtain the simple resul t 10 i 1:L i 15 i 19 

6k 

(V = constant, p) = F qi (p) \ Vi (q)= v (17) 

^Bulk 

for the Knudsen effect on the nanostructure thermal transport. This follows because the 
value of 7] is then independent of q and of the mode. For the case of wires and multilayered 
structures it is possible to obtain analytical results for the solution (|17|). 



III. INTERFACE-LIMITED NANOSTRUCTURES THERMAL TRANSPORT 

The study of the thermal transport in wires and layered structures requires a careful 
analysis 37 of the phonon scattering at the interfaces of the specific systems. From this 
analysis we can complete the evaluation of the Knudsen transport suppression (|17|) for a 
free-standing wire, for a individual layer with possibly different nature of scattering at the 
top and bottom interfaces, and for a multilayered structure. Details of the derivation is 
published elsewher e) 50 ' 51 and here we just summarize the approach and the analytical results 
obtained.— 
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Ill A. Thermal conductivity in wires 



For a free standing wire (Figure HJ we assume a thermal gradient along the wire-axis 
V. For simplicity we let the x-direction coincide with the perpendicular phonon momen- 
tum, i.e., q± = q x , and we assume the same specular scattering ratio p along the whole 
boundary. Under these conditions we solve the 1-BTE (Eq. (fTTjl ) for the change 6N W1IC 
subject to the boundary conditions and (jlOj) with f2 + = a/ (d/2) 2 — y 2 for x > and 
= —a/ (d/2) 2 — y 2 for x < 0. We analyze the boundary effect following Ref. |37|, adapt- 
j n gi5j50 description originally developed for electron transport in thin films. 

Using the definition ()12j) we formally determine the wire suppression ratio 



(18) 



^ire, / , = (1 x exp(-y/(d/2) 2 - y 2 / {r q v q (q x / q) }) p -, /{Wfe/g)} 

9 ^ 1 - pexp(-2 V(d/2)2 _ tf/favfa/q)}) 

and evaluate the average amplitude suppression^* 4 ^ 

K he (l±J^P) = -( 1 -p) 2 J^jP j ~ 1 / ^v/T^e^/^/ 9 '. (19) 

A polar coordinate system, where q x = q sin 9, is then used to express the weighted Knudsen 
suppression factor (JTfij) at mode i, 

F 7 C (P)\ Vl(g) = (y) j^^l^ 6 cos 2 dHJ c (q x /q = sind; p) 

oo 

= (i- p ) 2 J2jp } ~ 1 m(jv), (20) 



in terms of a Meijer's G special function^ and a polynomium: 

2 I 

(21) 



M(C) = ? ^2 



1,3 

1 3 _1 1 

2 ' 2 ' 2 ' 




For purely diffusive scattering (p —>■ 0) only the first term in the sum of Eq. (|20j) contributes 
to the lattice thermal conductivity and therefore the ratio « W ire/ ^Buik simplifies to 



K wire / „n 1 6 I 7T 3 ^ 2 ^2 / 3 

^Bulk 7? 4 



2' 2' 2' 

Similarly, if instead we assume pure specular scattering along the whole wire-boundary then 
F mre (p — > 1)| — > so that the bulk transport description is correctly recovered. 
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Our analytical results (|20I)- ([2"T]) and/or (}2*2*j) are in exact agreement with the numerical 
results presented in ReiU Our analytical 1-BTE result differs quantitatively but not 
qualitatively from the nonselfconsistent n wive determinatio n 18 ^ 9 that can be extracted from 
analytical calculations of how boundary scattering shortens the average phonon mean-free 
path.— We stress that our results (J2TH) - (}2"T]) and/or (j22j) completes the previously published 
investigation^ by a special function evaluation for a significant speed up in computation 
efficiency We also stress that care must be taken in the numerical evaluation of our special 
function (j2"Tj) . or the integral in Eq. ()20|) for large arguments (( > 10) when using, for 
example, the Mathematica program. 



Ill B. Thermal conductivity in an individual layer 



We consider a layer of thickness d with its boundaries at z = ±d/2 and a thermal gradient 
lying in the xy-plane. In this geometry we have q± = q z . In the general case, p + ^ p-, 
we must analyze the phonon scattering at both interfaces in order to evaluate the induced 
change, <5iVg ayer . Generalizing previous result a 10 i 1:L i 15 i 37 i 38 and introducing 77 = c?i ayei ./£ mfp and 
£ = rj/(2\q z \/q), we obtain for general the result^^i 



^ <i%Z>o(p+, P-) ^ 



V ^Sr<o(p+> p-) j 



1 



^ (p_ - l)e € +p_(p^ 



P-P+e 



-2* 



l)e 



-A 



\ (P+ - !)e € + p+(p- - l)e 5 y 



(23) 



for the amplitudes of the suppression caused by scattering at the upper and lower interfaces. 

For our single-layer and multilayer calculations we will generally assume that we have a 
constant £ m f p and can hence express the Knudsen transport suppression 



K Bulk 



) V— ) I r/=constant j 



(24) 



generalizations to a mode and frequency-dependent phonon mean-free path are as straight- 
forward as in the wire case. We calculate again the cross- section average of the distributions 
change, TC l g yer (q z /q;p + ,p-), and determine the phonon suppression factor as an average over 
the direction of momentum q 



F laycr (^ 



,P- 



1 77=constant 



jf d6 sin 3 9H l q aycr = cos 9; p+, p_J , (25) 
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noting that the sign of cos(#) = q±/q = q z /q uniquely determines whether the emerging 
phonons have scattered off the top or bottom material boundary 5 - Introducing finally 

P = P+ 2 V ~ and P = VP+P-- ( 26 ) 
and performing the ^-integration we obtain the analytical resul t 47 i 50 i 51 

F layer fep)|^ constant = l|(l-p)-4^g(i,r 7) p- p)J, (27) 

given in terms of a modified exponential error function 

Q{j, 7], p, p) = {p 2 - 2p+l)f- 1 D(jr ] ) + {p-p)(p+l) 2 p 2{i ' 1) D{2jr ] ) 7 (28) 

where 

y n \A 12 24 24 J 
At p + = p- the present more general result ()27j) - (j29j) correctly simplifies to 

^(v, P ) = Ji-^i-^W- 1 ^)}, (so) 

«Buik 8?7 ^ ^ j 

as obtained in previous der ivat ions . 10 > 1 1 > 15 > 32 

Figure H] shows the reduction of the lattice thermal conductivity in pure semiconductor 
layers. The top panel documents the ratio Kiayer/^Buik as a function of 77 = rfi ay er/^mfp for 
p+ = P- = P- As in the wire case, the lattice thermal conductivity naturally becomes higher 
to recover the bulk behavior when the probability for specular scattering increases. The 
lower panel shows /ti ayC r/ K Buik for a fixed layer-thickness <ii ayer = 0.1£ m f p and for varying 
specular scattering probabilities, p + and p_. The figure documents that the thermal con- 
ductivity remains suppressed below the bulk transport value even when the phonons are 
purely specularly scattered from one of the two interfaces. 



Ill C. Thermal transport in multilayered structures 

To estimate the thermal transport in heterostructures we make the approximation to 
ignore all phonon transmission between the layers. There exists exact results^ 3 - but our focus 

12 



is on providing simple estimates and extracting consequences of the size-limited thermal 
transport. We consider structures with layers j — 1, 2 ... of individual thicknesses dj and 
with probabilities for specular (rather than diffusive) scattering at the top (bottom) 

layer boundary. Our approximation then consists in approximating the effective (linear- 
response) conductivity as 

/'••/'' K' , Q1 s 

where Kj(dj,p 3 + ,p'_) denotes the thermal conductivity of layer j. We stress that the 
individual-layer thermal conductivity typically is subject to a significant Knudsen-flow sup- 
pression and must be calculated for the layer-specific nature of the boundary scattering 
an d f° r the specific individual-layer thickness dj, as also indicated by the argument 
to the individual-layer conductivities Kj. 



IV. RESULTS AND DISCUSSION 



The Knudsen flow effect has direct consequences for the semiconductor nanostructure 
thermal transport. Below we test the accuracy of our approximations and detail the size 
effects on the temperature variation of K wire and Ki aye r- 



IV A. Thermal transport in Si and SiC nanostructures 

Figure 03 tests the accuracy of our approximative determination of the Knudsen flow 
effects by comparing the measured 3 and calculated temperature variation of the lattice 
thermal conductivity in a 115 nm thick silicon wire. The solid circles represent experimental 
data from Ref. y| while the solid (dashed) curve shows our theoretical results (J2U|) - (flHj) . 

for the Knudsen flow calculated for pure diffusive scattering (p — 0) and assuming 
a frequency and polarization-independent bulk-phonon mean-free path £ m f p (a frequency- 
dependent relaxation time^ r q = AiT^joj 1 using ((7j) — (jHJ) to calculate the function A(T)). 
The insert compares the predicted Si wire thermal conductivity with measurements^ of bulk 
Si on logarithmic scales, emphasizing the significant effects of nanostructuring. 

The figure reveals some discrepancies which we assign to the simple approximation we 
use for the phonon dispersion.— The figure also shows that the more elaborate model with a 
frequency dependent r and £ m f p yields the more accurate description. However, since these 
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two models give similar magnitude and qualitative behavior of K w i re (with good agreement 
to the measured values) we find that the analytical result (}2"2")) serves as a simple, analytical 
estimate of the wire conductivity K w i re - 

Figure El compares the corresponding thermal- conductivity variation with temperature 
of bulk SiC (dashed-dotted curve) and of a 115 nm thick SiC wire (dashed curve) and 
layer (solid curve) structure. For both nanostructures we assume pure diffusive scattering 
(p = 0) at all interfaces and a frequency and mode-independent bulk-phonon mean-free path 
described by Eq. (JSJ). 

Together the insert of Fig. Eland Fig. distress the two features which are important when 
studying the thermal transport in nanostructure (Si and SiC) systems:—^ (i) the maximum 
bulk-transport value is about two orders of magnitude higher than the corresponding peak 
values for the nanostructures, i.e., «BuikC0 ~ 100 K iayer,wire(^)' an d (h) the maximum val- 
ues of the lattice thermal conductivity in the wires and layers are shifted towards higher 
temperatures from the peak of the bulk conductivities. 



IV B. Thermal transport in doped materials 

Figure [7| documents that the Knudsen suppression in the lattice thermal conductivity also 
affects the case of doped of Si and SiC nanostructures. The figure compares the conductivity 
of a Si and SiC layer at two different temperatures, 77 K and 300 K for pure materials 
(solid curves) and in the case of a very high dopant concentration ~ 10 20 cm -3 (dashed 
curves) . We assume here that the phonons are scattered purely diffusively at both interfaces 
ip+ — P- — 0) and use the analytical result P7 J) -(P§ )) for an assumed constant £ m f p (T). We 



use the measurements of Ref. |5y with the Si layer p-doped to be 1.7 • 10 20 cm -3 and the 
SiC layer n-doped 5 • 10 19 cm -3 . The phonon mean-free path in the doped layers is thus 
estimated at room temperature to be 49 nm and 138 nm for SiC and Si, respectively and at 
T = 77 K to be 420 nm and 830 nm, respectively. 



IV C. Breakdown of Fourier transport analysis 

Figure |H1 documents that traditional thermal transport calculations fail to accurately de- 
scribe the phonon transport in low-dimensional systems at and below room temperature. 
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This is illustrated for a Si02/SiC/Si02 heterostructure for which we have calculated the 
effective multilayer thermal conductivity (J3T|) for varying thickness of the SiC layer, (fee- 
At the SiC/Si02 interface we assume pure diffusive phonon scattering (p = 0). The solid 
curve represents the effective lattice thermal conductivity obtained by the present Knudsen 
approach (described in Sec. Ill B) in which the thermal-conductivity suppression in each 
individual layer is considered. In contrast, the dashed curve was obtained by replacing the 
suppressed value of the layer conductivity in Eq. (p!T]) by the corresponding bulk transport 
conductivity, as is typically done in many finite element calculations. This Fourier analysis 
yields an effective lattice thermal conductivity which is too large. Traditional Fourier anal- 
ysis 39 of the lattice thermal transport is inadequate for some nanostructured semiconductor 
systems in the low-temperature regime and the Knudsen-flow suppression affects even the 
room-temperature behavior. 

IV D. SiC oxidation 

Figure M illustrates the use of the Knudsen-flow thermal transport as a possible gauge of 
the SiC oxidation process. Metal-oxide-semiconductor transistor gates are usually formed 
by oxidizing the Si to form a thin oxide filin-S located on top of the pure transport channel. 
This has been the focus of many experimental and theoretical investigations .^ 1 ^ 2 ^ 3 ^ 41 ^ 1 ^ 
Many of the experimental probes require that one extracts the sample and investigate it 
in, for example, a transmission electron microscope whereas qualitative in-situ probes could 
be more useful for some applications. We note that the SiC/Si0 2 interface is rough but of 
relative uniform thickness.— The situation is similar to the case of adsorption on metal films 
where the (electrical) conductive is proven to be sensitiv o 67 i 68 i 69 i 70 i 71 i 72 to the film thickness 
and surface nature. We therefore argue 4 ^ that measurements of the thermal Knudsen flow 
could serve as such an in-situ probe of the oxidation progress .^ 1 ^ 2 ^ 3 ^ 41 ^ 1 ^ 

The figure shows the variation of the effective lattice thermal conductivity, « e g, with the 
oxide-thickness, d Si02 , calculated by Eq. (pTTj). We use the value 7 - 3 - 1.4 W m _1 K _1 for the 
bulk thermal conductivity of SiC>2 to calculate the effective lattice thermal conductivity in 
two different heterostructures. In addition, we assume a uniform oxide-thickness during the 
whole oxidation process. The dotted curve corresponds to a single SiC layer on which a thin 
uniform Si02 film grows. For simplicity we assume that p = at the Si02/SiC interface and 
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that the thickness <isi0 2 + ^Sic = 200 nm is kept fixed during the whole oxidation process 
that we model.— The solid (dashed) curve correspond to results for a doped-SiC/SiC/Si02 
heterostructure evaluated under the assumption that p — 1 (p = 0) characterizes the scatter- 
ing at the doped-SiC/SiC interface; we also assume that the underlying SiC layer is doped 
at a concentration of 5 • 10 19 cm -3 . Our results document that the effective lattice thermal 
conductivity is sensitive to the oxide thickness and could be used as a qualitative in-situ 
probe of the SiC oxidation progress.— 

V. CONCLUSIONS 

We have theoretically investigated the in-plane thermal transport in nanostructured semi- 
conductor wires, layers and multilayers, in the framework of the 1-BTE. The calculations 
include the important boundary scattering which specifies the phonon Knudsen flow and 
plays a dominant role in determining the thermal transport in semiconductor nanostruc- 
tures. We obtain analytical results for the nanostructure thermal conductivity under the 
assumption of a frequency independent (constant) mean-free path £ mfp . Our results further 
simplifies the thermal-transport calculations under more accurate assumptions for £ m f p . We 
find that the lattice thermal conductivity is significantly suppressed from the bulk value at 
room temperature and below in both thin-wire and thin-layer nanostructures. 

We have tested our approximative description of the phonon Knudsen flow effect by 
comparison to measurements of the thermal conductivity for thin Si wires and find that our 
analytical result for a constant £ m f p provides a fair accuracy in a qualitative description. 
We further show that the traditional Fourier analysis^ 9 - (typically used in finite-element 
calculations) predict an effective lattice thermal conductivity which is inaccurate for some 
nanostructured semiconductor systems especially in the low-temperature regime. Finally we 
have illustrated the use of the Knudsen-suppression thermal transport as a possible gauge 
of the SiC oxidation process. 
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TABLES 



TABLE I: Model parameters of the phonon description in SiC, Si, GaAs, and Ge. The measured 



longitudinal and transverse sound velocities, q and q, in SiC are taken from Ref. 
Si, GaAs, and Ge we use Ref. 



53 while for 



52l . The Debye temperatures, G; and 0t, are defined as 0j 



Ref. 



a X )/fc_B and have in all four cases been calculated from the phonon dispersion curves in 



Material 


Q 






e t 


cxiuf(q max ) 


a^ 3 (g max ) 


SiC 


13300 


7250 


880 


380 


0.72 


1.23 


Si 


8480 


5870 


570 


160 


0.43 


1.91 


Ge 


4960 


3570 


360 


120 


0.20 


1.51 


GaAs 


4730 


3340 


270 


100 


0.55 


1.95 
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FIG. 1: Typical semiconductor nanostructures that are subject to a Knudsen-flow effect (Refs. 
and [3^ ) with a suppression of the thermal conductivity. Phonon transport can be effectively 
confined to device channels (light gray) by the low conductivity of the processed surrounding. 
Moreover this in-channel thermal transport typically suffers a significant suppression by interface 
scattering. The upper panel illustrates two possible realizations of a wire which could be embedded 
in an oxide or free-standing. Phonons experience specular (diffusive) scattering with the probability 
p (1—p) at the wire boundaries. The lower panel shows a corresponding realization of a multilayered 
nanostructure in which, for example, a transport channel is located between two (in general) 
different materials (for example, a top oxide and a heavily doped substrate). Phonons scatter 
specularly from the upper (lower) interface with probability P+i-) and diffusively with probability 
1 — P+(~) to produce a Knudsen flow. 
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FIG. 2: Summary and origin of predicted suppression in the thermal conductivity of nanostrucured 

semiconductor structures. The left panels show the thermal conductivity of layer (solid curves) and 

wire (dashed curves) structures calculated as a function of the layer thickness and wire diameter 

at room temperatures for three typical semiconductor materials. The dotted lines identifies the 

experimentally observed value of the bulk-material thermal conductivities from which we also 

extract the listed values of the phonon mean free paths £ m f p . A corresponding graph for Ge 

gives almost similar values to that of GaAs. The right top and right middle panels compare the 

predicted temperature variation of the layer and wire thermal conductivity in SiC and Si for a 

layer thicknesses di^er/wire diameter d w ; re for which there now exist experimental realizations as 

Si nanowires (Ref. 3). Finally, the lower right panel shows the calculated temperature variation in 

the almost micron-scale phonon mean free paths in bulk SiC, Si, and GaAs. The predicted phonon 

mean free path exceeds the typical separation between material interfaces in nanostructures and 

explains the device importance of the phonon Knudsen flow effect. 
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FIG. 3: Schematics and details of the Knudsen-flow effect as expressed in the momentum and 
spatial variation in the thermal-gradient induced changes in the phonon distribution. The figure 
shows (for a given magnitude q of the phonon momentum q) the ratio of the distribution change 
8N q = <5iV™ ire (arising in a wire with purely diffusive interface scattering) to the corresponding 
bulk-distribution changes 5iV? evaluated for different ratios of cf w i re /^ m fp(<7)- The top panel 
shows this distribution ratio (averaged over the wire cross section) as a function of the angel 9 
between the phonon momentum and the wire axis/thermal-gradient direction/transport direction 
z. The bottom panel shows the corresponding variation of the distribution ratio (averaged over 9) 
with the coordinate x (at y = 0) along the wire diameter which is assumed to coincide with the 
direction of the perpendicular momentum component q±. 



24 




FIG. 4: Relative suppression of the thermal conductivity in a general individual semiconductor 
layer evaluated for a the constant mean free path 4nfp- The upper panel shows the thermal conduc- 
tivity ratios, ftiaycr/^Buik, as a function of r\ for specular scattering probabilities p + = p_ = p. The 
thermal conductivity naturally increases to recover the bulk-like behavior recovers when the prob- 
ability for specular scattering increases. The lower panel shows the ratio Ki ayer /KBuik as a function 
of p + for different p_ evaluated at a given ratio 77 = diayer/Anfp = 0.1. The panel documents that 
the layer thermal conductivity remains significantly suppressed from the bulk transport value even 
when the phonons scatter purely specularly (p± = 1) from one of the interfaces. 
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FIG. 5: Test of accuracy of our approximative phonon Knudsen flow description. The figure shows 
the temperature variation of the lattice thermal conductivity in a 115 nm thick silicon wire. The 
solid circles represent experimental values from Ref. |^. The dashed (solid) curve shows our theory 
estimate for pure diffusive scattering, p = 0, and calculated assuming a frequency-dependent 
relaxation timei^ r q = A(T)/uj^ (a frequency-/polarization-independent bulk-phonon mean free 
path £ m fp(T) specified by bulk- measurements.) The inset compares the temperature-variation on 
logarithmic scales of bulk-Si thermal transport (dashed-dotted curve) and the evaluated Si-wire 
conductivity (solid curve) calculated for a frequency- independent £ m f p (T) . 
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FIG. 6: Thermal conductivity suppression in 115 nm thick SiC wires (dashed curve) and lay- 
ers (solid curve) compared to the bulk-SiC thermal conductivity measurements (dashed-dotted 
curve) fSfi The peak values for the wire and layer transport case are about two orders of magnitude 
smaller than the bulk-transport peak and are shifted to much higher temperatures. 
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FIG. 7: Reduction of the lattice thermal conductivity in doped and undoped layers of Si (top 
panels) and SiC (bottom panels). In all four panels the lattice thermal conductivity, Ki aycr , is 
plotted as a function of the layer-thickness, cZi ayer , for p + = p_ = at T = 300 K (left panels) and 
T = 77 K (right panels). The top panels show /q aye r in pure (solid curves) and in p-doped (dashed 
curves) Si layers. Similarly, the bottom panels compare the lattice thermal conductivity in pure 
SiC layers (solid curves) to that of n-doped layers (dashed curves). 
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FIG. 8: Comparison of traditional thermal-transport calculations 3 ^ (dashed curve) and present 
Knudsen-flow results (solid curve), Eqs. (|3T|) and H24 |) -(|29 |l . The panel shows the effective lattice 
thermal conductivity as a function of the layer thickness in a SiC layer that is located between 
two 50 nm thick SiG"2 layers. The traditional thermal-transport calculations overestimates the 
thermal conductivity significantly in low-dimensional semiconductor structures at room tempera- 
tures and below. The figure therefore suggests a breakdown of the Fourier approach based on bulk 
conductivities as in typically used in finite-element calculations. 
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FIG. 9: Possible probe of the SiC oxidation process as measured by the variation in the effective 
lattice thermal conductivity with the oxide-thickness, dsio 2 - The dotted curve correspond to a 
single SiC layer on which a thin uniform SiO"2 film grows. The solid (dashed) curve shows results 
for a doped-SiC/SiC/Si02 heterostructure evaluated under the assumption that p = 1 (p = 0) 
characterize the scattering at doped-SiC/SiC interface and p = characterize the SiC^/SiC inter- 
face. We assume that the underlying SiC layer is doped at a concentration of 5 • 10 19 cm -3 . For 
simplicity we also assume that the thickness dsio 2 + ^SiC = 200 nm is kept fixed during the whole 
oxidation process that we modeU*. 
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